Interplay between phonon and impurity scattering in 2D hole transport 
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We investigate temperature dependent transport properties of two-dimensional p-GaAs systems 
taking into account both hole-phonon and hole-impurity scattering effects. By analyzing the hole 
mobility data of p-GaAs in the temperature range 10 K< T <100 K, we estimate the value of the 
appropriate deformation potential for hole-phonon coupling. Due to the interplay between hole- 
phonon and hole-impurity scattering the calculated temperature-dependent resistivity shows inter- 
esting nonmonotonic behavior. In particular, we find that there is a temperature range (typically 2 
K< T <10 K) in which the calculated resistivity becomes independent of temperature due to a sub- 
tle cancellation between the temperature dependent resistive scattering contributions arising from 
impurities and phonons. This resistivity saturation regime appears at low carrier densities when 
the increasing resistivity due to phonon scattering compensates for the decreasing resistivity due 
to the nondegeneracy effect. This temperature-independent flat resistivity regime is experimentally 
accessible and may have already been observed in a recent experiment. 



I. BACKGROUND 



The room-temperature resistivity of a metal (as well 
as most metallic electronic materials, e.g., doped semi- 
conductors) is mainly limited by the electron-phonon 
interaction [1| i.e., by phonon scattering, with a no- 
table exception being doped graphene with its very weak 
electron-phonon coupling [2]. The electron-phonon scat- 
tering contribution to the resistivity falls off strongly 
at very low temperatures (T < Tbg) in the so-called 
Bloch-Griineissen (BG) regime [1], due to the exponen- 
tial suppression of the bosonic thermal occupancy of the 
phonons. At low temperatures, therefore, the electronic 
conductivity of metallic systems is invariably limited by 
disorder, i.e., by electron-impurity scattering, which gives 
rise to the zero-temperature residual resistivity of metals 
at low temperatures where the phonon scattering contri- 
bution has vanished. In general, the impurity scattering 
contribution to the metallic resistivity is temperature in- 
dependent, at least in three-dimensional (3D) metals, be- 
cause the impurities are quenched, and the temperature 
scale is, therefore, the Fermi temperature (Tp ~ 10^ K in 
3D metals), the impurity scattering contribution to the 
resistivity is essentially temperature-independent in the 
— 300 K regime. This, however, is not true if Tp is low as 
it could be in 2D semiconductor-based systems with their 
tunable carrier density where at low densities Tp could be 
just a few K. In fact, for 2D GaAs hole systems (2DHS), 
the situation with Tp < 1 K can easily be reached for a 2D 
hole density ~ 10^^ cm~^ [31,05. This situation, which has 
no known analog in 3D metallic systems, leads to very 
strong experimentally observed temperature-dependent 
2DHS resistivity in the T < 1 K temperature range aris- 
ing entirely from electron-impurity scattering since the 
phonon scattering is completely thermally quenched at 
such low temperatures. As an aside, we note that the 
temperature dependence of 2D electronic resistivity aris- 
ing from impurity scattering in 2D GaAs based electron 



systems is rather weak due to the much lower effective 
mass of 2D electrons {rrie ~ O.OTttiq, where mo is the 
vacuum electron mass), leading to much higher Tp, com- 
pared with 2D holes (mh ~ 0.4mo) in GaAs-based 2D 
heterostructures. 

The goal of the current work is to explore the interplay 
between impurity scattering and phonon scattering in 
the temperature-dependent resistivity of 2DHS in GaAs- 
based 2D systems. At higher temperature (> 100 K), 
the GaAs carrier resistivity is completely dominated by 
longitudinal optical (LO) phonon scattering, which has 
been extensively studied [5] and is not a subject matter of 
interest here since the resistivity limited by LO-phonon 
scattering manifests the strong exponential temperature 
dependence ^ g-^c^Lo/fceT ^[j^i^ huji^o ~ 36 meV. 

Our interest is in the interplay between acoustic 
phonon scattering and impurity scattering in the low to 
intermediate temperature regime (T ~ Tbg — ^f) where 
both impurity scattering and acoustic phonon scattering 
contributions to resistivity would show nontrivial tem- 
perature dependence. In particular, we are interested in 
the question whether an interplay between the two scat- 
tering processes could lead to an approximately constant 
(i.e. temperature independent) resistivity over some in- 
termediate temperature range (1 K < T < 40 K). One 
of our motivation comes from a recent experiment p 
which discovered such an intermediate-temperature "re- 
sistivity saturation" phenomenon in a 2D GaAs-based 
electron system in the presence of a parallel magnetic 
field (which presumably serves to enhance the electron 
effective mass due to the magneto-orbital coupling, thus 
reducing the Fermi temperature of the electron system 
[7]). A second motivation of our work is estimating the 
deformation potential coupling strength for hole-acoustic 
phonon scattering in 2DHS in GaAs. It turns out that 
the electron-phonon deformation coupling is not known 
in GaAs, and a quantitative comparison between our the- 
oretical results and experimental transport data could 



lead to an accurate estimation of the deformation poten- 
tial coupling in the GaAs-based 2DHS. We mention in 
this context that the accurate evaluation of the electron 
deformation potential coupling in 2D GaAs systems is 
also based on a quantitative comparison of the experi- 
mental and theoretical transport data [5|, [8| . 

The basic physics we are interested in (see Fig. 1) is 
a situation where the acoustic phonon contribution to 
the resistivity increases linearly with increasing tempera- 
ture (T > Tbg), but the impurity contribution decreases 
with increasing temperature (T > Tp) as happens in a 
nondegenerate classical system where increasing temper- 
ature must necessarily increase the conductivity since the 
electrons are classically moving "faster" . (We emphasize 
that such a situation is physically impossible in 3D met- 
als since T > 10^ K would be required, but is routinely 
achieved in 2D semiconductor-based systems where very 
low carrier density could lead to very low values of Tp.) 
We explicitly consider a situation with T <C Tlo, where 
Tlo ~ 100 K in GaAs where LO-phonon start contribut- 
ing to the resistivity in a substantial manner. Now we 
ask the question whether it is possible for the increasing 
temperature dependent contribution to the hole resistiv- 
ity arising from enhanced phonon scattering at higher 
temperatures could be approximately canceled over a fi- 
nite temperature range by the decreasing temperature 
dependent contribution arising from impurity scattering 
with increasing temperature due to the quantum-classical 
crossover, i.e., the nondegeneracy, effects. Some early 
theoretical work |9] indicated that the interplay between 
phonon scattering and nondegeneracy may indeed lead to 
a partial cancellation of different temperature-dependent 
contributions, producing an interesting nonmonotonicity 
(see Fig. 1) in the resistivity as a function of temperature 
in the 1 — 10 K range for low-density, high- mobility 2DHS 
in GaAs. In the current work, we look into this issue at 
great depth in view of the recent experimental work. In 
addition, we obtain the appropriate 2DHS deformation 
potential coupling constant by comparing our theoretical 
results to existing 2DIIS hole transport data. 



II. INTRODUCTION 

The temperature dependent transport in 2D systems 
has been a subject of intense activity since the obser- 
vation of an apparent electronic metal-insulator transi- 
tion (MIT), which represents the experimental observa- 
tion of a transition from an apparent metallic behavior 
(i.e., dp/dT > 0, where p is the 2D resistivity) to an insu- 
lating behaviors (i.e. dp/dT < 0) as the carrier density is 
reduced. The remarkable observation of the anomalous 
metallic temperature dependence of the resistivity is ob- 
served mostly in high-mobility low-density 2D semicon- 
ductor systems [lO, lll|. The low temperature anoma- 
lous metallic behavior discovered in 2D semiconductor 
systems arises from the physical mechanism of strong 
temperature dependent screening of charged impurity 



scattering [lO]. At low temperatures (T < Tp, where 
Tp = Ep/kB is the Fermi temperature with the Fermi 
energy E-p) the main scattering mechanism in resistiv- 
ity is due to impurity disorder from unintentional back- 
ground charged impurities and/or intentional dopants in 
the modulation-doping. The resistivity Pimp(^) limited 
by the charged impurities increases linearly with temper- 
ature at lower temperatures (T < Tp) due to screening 
effects. This is a direct manifestation of the weakening 
of the screened charged disorder with increasing temper- 
ature [10, [l2| or equivalently an electron-electron inter- 
action effect in the so-called ballistic regime [13|. For 
T ^ Ty^ p{T) decreases as Tp/T due to nondegeneracy 
effects and the quantum-classical crossover occurs at the 
intermediate temperature regime around T ^ T^. [l2| 

The temperature dependent resistivity [5| limited by 
acoustic phonons pph(T) undergoes a smooth transition 
from a linear-T dependence at high (T > Tbg) temper- 
atures to a weak high power T^ dependence with a > 3 
as the temperature is reduced below /cbTbg = 2/cFVph, 
where k^ is the Fermi wave vector of the 2D hole system 
and v^\^ is the longitudinal or transversal sound veloc- 
ity. The characteristic temperature Tbg is referred to as 
the Bloch-Griineisen (BG) temperature. Note that Tbg 
is much smaller than the Debye temperature since the 
inverse lattice constant greatly exceeds /cp. In the BG 
regime (T < Tbg) the scattering rate is strongly reduced 
by the thermal occupation factors because the phonon 
population decreases exponentially and the phonon emis- 
sion is prohibited by the sharp Fermi distribution, which 
gives rise to high power law behavior (a > 5 with screen- 
ing effects, but a > 3 without screening) in the tem- 
perature dependent resistivity. Thus, the phonon con- 
tribution to the resistivity is negligible compared to the 
charged impurities and the phonon contribution to the 
resistivity shows very weak temperature dependence for 
T < Tbg- For temperatures T > Tbg, since the electron- 
phonon scattering becomes proportional to the square of 
the oscillation amplitude of ions, the pph depends linearly 
on the temperature. We note that at low carrier density 
where /cp is small, Tbg can be very low. 

In this paper we investigate the temperature depen- 
dent transport properties of p-type GaAs-based 2DHSs 
for the temperatures T < 100 K by taking into ac- 
count both hole-phonon and hole-impurity scatterings. 
In p-GaAs, holes interact with acoustic phonons through 
a short-range deformation potential as well as through 
a long range electrostatic potential resulting from the 
piezoelectric effect. The precise value of the deformation- 
potential coupling constant D is very important to un- 
derstand 2D p-GaAs carrier transport properties. For 
example, the mobility is strongly related to the defor- 
mation potential, p ~ D~^ . However, the precise value 
of D for p-GaAs has not been available, so the value 
of the n-GaAs deformation potential (T) = 12 eV Q) 
is used uncritically. Thus, the current investigation of 
the transport properties of p-GaAs systems is motivated 
by getting the exact value of the deformation potential 



constant D in p-GaAs. We find that the fitted value of 
deformation-potential coupling could change by as much 
as 60% (i.e. D = 7.6 — 12.7 eV) depending on the value of 
the depletion density ndepi- When we assume ndepi = 0, 
we obtain D = 12.7 eV as the most suitable value for 
the p-GaAs acoustic phonon deformation coupHng con- 
stant by fitting the available experimental data [14]. The 
value of D = 12.7 eV is larger than the generally ac- 
cepted value in bulk GaAs {D = 7 eV) [l5| but compa- 
rable to the value of the n-GaAs (12 - 14 eV) [i,[ii,[l3- 
However, due to the lack of knowledge about ndepi we 
have the uncertainty in the value of deformation poten- 
tial coupling. To precisely determine the value of D for 
holes it is required a careful measurement of p{T) over the 
T = 1 — 50 K range in a high mobility and high density 
(n > 5 X 10^^ cm~^) sample. Note that our finding of the 
apparent dependence of the GaAs 2D hole deformation 
potential coupling on the background (and generally un- 
known) depletion charge density in the 2DHS is obviously 
not a real effect since the hole-phonon coupling strength 
cannot depend on the depletion charge density. The ap- 
parent dependence we find arises from the fact that the 
calculated resistivity depends strongly on the depletion 
density whose precise value is unknown. 

Other motivation of our work is to investigate the non- 
trivial nonmonotonic behavior of the temperature depen- 
dent resistivity observed experimentally. In the pres- 
ence of the hole-phonon and hole-impurity scattering a 
nonmonotonicity arises from a competition among three 
mechanisms [j, 0, [isi]: screening which is particularly 
important for T < T^^ nondegeneracy and the associ- 
ated quantum-classical crossover for T ^ T^^ and phonon 
scattering effect which becomes increasingly important 
for T > 5 — 10 K, depending on the carrier density. In 
Fig. [1] we show the schematic resistivity behavior of the 
p-GaAs in the presence of both charged impurity and 
phonon scatterings. At low temperatures (T < Ti ~ Tp), 
the scattering arising from charged impurities dominates 
and the resistivity increases as temperature increases due 
to the screening effects. For T > T2 ~5 — lOK the scat- 
tering by acoustic phonons plays a major role and limits 
the carrier mobility of p-GaAs systems in this temper- 
ature range. Note that in general T2 :^ Tbg for p- 
GaAs systems. However, at intermediate temperatures 
(Ti < T < T2) the competition between acoustic phonon 
scattering and impurity scattering gives rise to a nontriv- 
ial transport behavior. We carefully study the non-trivial 
transport properties of p-GaAs in the intermediate tem- 
perature range (i.e., Ti < T < T2). Interestingly we 
find that the approximate temperature dependent p{T) 
in an intermediate temperature range can be constant, 
which arises from the approximate cancellation between 
the quantum-classical crossover and phonon scattering, 
and is quite general. When the phonon scattering dom- 
inates impurity scattering for T > T2 the temperature 
dependence of hole mobility enables one to extract infor- 
mation on the electron-phonon scattering from mobility 
measurements. In this temperature range the valence 



FIG. 1. (Color online) (a) Schematic resistivity behavior in 
the presence of both hole- impurity scattering and acoustic 
phonon scattering as a function of temperature. Typically 
Ti ^ Tf and T2 ^ 5 — 10 K depending on the carrier density. 



band deformation potential can be determined by fitting 
theoretical calculations to the existing carrier mobility 
data. In this paper we extract the value of the deforma- 
tion potential by fitting the experimental mobility data, 
but uncertainty arises from the unknown values of the de- 
pletion charge density and background impurity charge 
density. 

To investigate the temperature dependent transport 
properties for the p-GaAs 2DHS we use the finite temper- 
ature Boltzmann transport theory considering scatter- 
ings by charged random impurity centers and by electron- 
acoustic phonons with finite temperature and finite wave 
vector screening through random phase approximation 
(RPA) [10]. We also include the finite size confinement 
effects (i.e., we take into account the extent of the 2D 
system in the third dimension and do not assume it to 
be a zero- width 2D layer). The effect we neglect in our 
theory is the inelastic optical phonon scattering. Polar 
carrier scattering by LO phonons is important in GaAs 
only at relatively high temperatures (T > 100 K), be- 
coming dominant at room temperatures. Due to the 
rather high energy of the GaAs optical phonons (~ 36 
meV), LO phonon scattering is completely suppressed in 
the temperature regime (T < 100 K) of interest to us in 
this work. Resistive scattering by optical phonons in 2D 
GaAs system has been considered in the literature [5[- 
We note that the other scattering mechanisms (e.g., in- 
terface roughness scattering and alloy disorder scattering 
in Gaa^Ali_ajAs, etc.) are known to be much less quanti- 
tatively important [19|] than the mechanisms we are con- 
sidering (i.e., impurity scattering and acoustic phonons) 
in this work. 

The rest of the paper is organized as follows. In Sec. 
HI we present the general theory of the impurity and 
phonon scatterings, and discuss the power-law behavior 
of the resistivity in the low and high temperature limit. 
In Sec. IV we show our calculated resistivity taking into 



account hole-phonon and hole- impurity scatterings, and 
demonstrate that there is a temperature-independent re- 
gion due to the competition between the two scatterings. 
Finahy we conclude in Sec. V summarizing our results. 



III. THEORY 

To investigate the temperature dependent resistivity 
p(T) (or equivalently conductivity cr{T) = p{T)~^) of p- 
GaAs systems we start with the Drude-Boltzmann semi- 
classical formula for 2D transport. Due to the finite 
extent in the z-direction of the real 2D semiconduc- 
tor system we have to include a form factor depend- 
ing on the details of the 2D structure. In GaAs het- 
erostructures, the carriers are spatially confined at the 
2D interface and there is no longer translational invari- 
ance along the direction normal to the interface, desig- 
nated as the z direction. We assume that the confine- 
ment profile is described by the variational wavefunction 
^(x,^,z) = ^Q(^)e^(^-^+^^?/) [195, where 



^o{z) 



-b^z'^exp f "2^^ ) ' 



(1) 



and 6 is a variational parameter. For a triangular well, b 
is given by [20j 



b = 



/487rm*e2 



ridepl + ^^ 



(2) 



where m* is an effective mass, eo is the dielectric con- 
stant, ndepi is the depletion charge per unit area and n 
is the 2D carrier density. 

The density of states of 2DIIS with parabolic energy 
dispersion e{k) = l^ is given by 



D{s) 



Do {e>0), 
(e < 0), 



(3) 



The finite temperature Thomas-Fermi (TF) screening 
wavevector qTF{T) is defined by 



qTF{T) 



27re^ dn 



qtf 



1 — exp 



=1 



ksTj 



(7) 



where gxF = ^Do = ^i^- Note that 

r<ZTF (T^O), 



(8) 



A. Electron-phonon interactions 

The electron-longitudinal acoustic phonon interaction 
Hamiltonian for a small phonon wavevector Q is 



Hj,p{Q) = DQ-SR, 



(9) 



where D is the acoustic phonon deformation potential 
and SR is the displacement vector. Note that position 
operator in a simple harmonic oscillator with mass m and 
angular frequency uj is given by 



h 



2muj 



{a^a^), 



(10) 



where a and a^ is the annihilation and creation operators, 
respectively. Similarly, SR for Q can be expressed in 
terms of phonon annihilation and creation operators as 



5R{Q) 



2pm^Q 



eg (aq + a^g j , 



(11) 



where p^ is the mass density, ujq = viQ^ vi is the longi- 
tudinal sound velocity and eg is the phonon polarization 
unit vector. Thus the electron-phonon interaction by the 
deformation potential coupling can be expressed as 



where I^o = ^^ and ^g = 2 is the spin degeneracy 
factor. The carrier density n is given by 



Do 
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where f{s) = [e'^'-^-t' 



(4) 



1]~^ is the Fermi distribution 
function and P = 1/(/cbT). Alternatively, the chemical 
potential p{T) at a finite temperature T can be expressed 
as 



p{T) = kBTln [exp {E^/kBT) - 1] , 



(5) 



where ^f = tl/Dq. Note that Wvht^q piT) = Ey- From 
Eqs. (jl]) and dHj), 
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p{Q) (12) 



where p{Q) is the electron density operator. 

For the piezoelectric scattering in polar semiconduc- 
tors (i.e. GaAs) the scattering matrix elements are equiv- 
alent to the substitution of the deformation with the fol- 
lowing form [21[ 



D^^ 






(13) 



where /114 is the basic piezoelectric tensor component and 
A is a dimensionless anisotropy factor that depends on 
the direction of the phonon wave vector in the crystal 
lattice. We provide details of the parameter A in the 
following subsection. 



B. Boltzmann transport theory 

We calculate the temperature dependence of the hole 
resistivity by considering screened acoustic-phonon scat- 
tering. We include both deformation potential and piezo- 
electric coupling of the 2D holes to 3D acoustic phonons 
of GaAs. Details of the acoustic-phonon scattering the- 
ory are given in Ref. [55 . 

The transport relaxation time r{ek) at an energy Sk 
and a 2D wavevector k is given by 

-^= [ ^Wk,k, (1 - COS <p^^, ) \ ~ {^'"'^ , (14) 

T(£fc) J '^"'^ 



(27r)2 



1 - /(efc) 



where Wk,k' is the scattering probability between k and 
k^ states, and (j)kk' is the scattering angle between k and 
k' vectors. 

First, consider impurity scattering. Assume that im- 
purity charges are distributed randomly in a 2D plane 
located at (— (iimp)^ with a 2D impurity density nimp- 
Then the effective impurity potential for a 2D wavevec- 
tor q is given by 



/»oo 
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-qidimp-\-z) 



eo<7 



-qdimp 



b^q 



(15) 



From the Fermi's golden rule, Wk,k' for the impurity scat- 
tering has the following form: 



2^\^rnAQSl 



TTT-imp 



nimpS{ek -£k)^ (16) 



where q = k — k' . The dielectric function e{q^T) takes 
into account the screening effects of electron gas at a 
wavevector q and a temperature T. We will consider 
screened scattering within RPA approximation defined 

by 



eiq,T) = l 



qsiq,T) 



(17) 



where q^{q^T) is the temperature- and wave vector- 
dependent screening wavevector [22|. In the long wave- 
length limit (g = 0), qs{q,T) is given by Eq. ([7]). 

For electron-phonon scattering, Wk,k' for the electron- 
longitudinal acoustic phonon interaction has the follow- 
ing form [1]: 

^.,.^ = yy^^2(^;^A(..,.,)|/fe)| , (18) 
where Iz{qz) is the wavefunction overlap at qz defined by 
Iz{qz) = ldzeo{z)e^'^'. (19) 

Note that from Eq. ([I]), \Iz{qz)\^ = jE^f^' 



The factor C{q,qz) is the matrix element for acoustic 
phonon scattering. From Eq. (pT]) . for the deformation 
potential (DP), 



\CDp{q,qz 



,2 ^ D^^Q 
2p^vi ' 



(20) 



while for the piezoelectric potential (PE) [5], 

\c..M,,q.)\'=^^^^4^^^, (21) 



where Q = {q^qz)^ v\ {vt) is the longitudinal (transverse) 
sound velocity and 



9q'^q^ 



iW^cf 



(22) 



The factor A (£,£:'') is given by 

/\{e,e') = Nq5{e-e' ^fyjOq)^{Nq^l)5{e-e' -hwq), (23) 

where Nq = (^e^^^^ ~ 1)~^ is the phonon occupation num- 
ber. Note that the first and second terms in Eq. (|23|) 
correspond to absorption and emission of phonons, re- 
spectively. 

Finally, the total transport relaxation time is given by 

1 1 1 1 o 

, (24) 



Ttot(e) Tinip(£) TDp(e) TpE,l(£) TpE,t(£) 

where in the last term the degeneracy of the transverse 
modes has been taken into account. Then electrical con- 
ductivity in a 2DHS is given by 



2 [ dPk vl , .( df 



(25) 



where Vk is the mean- velocity at k. By inverting the 
conductivity we have the resistivity, i.e. p{T) = a~^{T). 



C. Quasi-elastic limit 

For a degenerate system (/cbT <C ^^f), all the scat- 
tering events at an energy s take place in a thin shell 
around the energy circle and the scattering can be con- 
sidered as quasi-elastic. For electron-phonon scatterings, 
the transport relaxation time in Eq. (p!4|) is given by 



27r f del) f dqz \C{q)^ \\'^ nt \f-\ a\ 

— Do I — / :^^^\I{qz)\ G((7)(l-cos(^), 



T{e) h V 27r7 27r e^{q) 

where q = 2/cf sin((/)/2) and G{q) is given by 

G{q) =/3j def{s) [(1 - f{e + hwg)) iV, 
+ (1 - /(£ - fuv,)) {N, + 1)] 

= 2ptUVgNg{Ng + 1). 



(26) 



(27) 



Thus for DP, Eq. (J26j) becomes ^ Q 
1 SD'^m'^bkBT r dcj) (1 
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while for PE, 
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where q = 9/32, Ct = 13/32 and 

1 + 6i^ + 13^2 + 2i^3 



/iH 
AH 



(1 + ^)6 

13 + 78i^ + 721^2 



(30) 



821^'=^ + 36i(;'^ + ^w^ 



13(1 



For simphcity, consider fx{q/b) = 1 case assuming b :^ 
1 in the extreme 2D hmit. Then, in the high temperature 
hmit, 



1 



3D'^m*bkBT 



TDP {T) 
1 



ocT, 



Wh^PmVf 

CA(e/ii4)^m* 4 T 



(31) 



(xT, 



rpE,A(T) 2n3p^^2 ^^^^ 

while in the low temperature limit, 
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where xtf = <7tf/(2/cf). Note that the resistivity is pro- 
portional to the inverse relaxation time, thus it follows 
the same power-law dependence in the high or low tem- 
perature limit. 

For the charged impurity with temperature-dependent 
screening wave vector the asymptotic low- and high- tem- 
perature behaviors of 2D resistivity for a S-ldiyei system 
are given by [23.] 



,(T < Tf) - po 



2xtf T 






,(r > Tf) ~ pi 



Tf 



1 + xtf Tf 

3^/^TXTF fTp^ 



3/2' 



(33) 
, (34) 



where po = p(T = 0), pi = (/i/e^)(nimp/n7rx|p), 
and C = 2.646 [xtf/(1 + ^tf)]^- At low temperatures 
(T < Tp) the resistivity increases linearly due to screen- 
ing (or electron-electron interaction) effects on the impu- 
rity scattering [10!,[l2|- To get the linear temperature de- 
pendent resistivity at low temperatures it is crucial to in- 
clude the temperature-dependent screening wave vector 



[1Q| |. At high temperatures (T > Tp) p{T) decreases in- 
verse linearly due to nondegeneracy effects. Thus, it is ex- 
pected that the resistivity has a maximum value and the 
quantum-classical crossover occurs at the intermediate 
temperature regime around Tp. When we consider both 
hole-phonon and hole-impurity scatterings the tempera- 
ture dependent resistivity becomes nontrivial due to the 
competition between these independent two scattering 
mechanisms. Since the resistivity limited by charged im- 
purities decreases at high temperatures phonon scatter- 
ing eventually takes over and p{T) increases with T again, 
which gives rise to nonmonotonicity in p(T). The non- 
monotonicity becomes pronounced in the systems with 
strong impurity scattering or at low carrier density. For 
weaker impurity scattering the phonon scattering domi- 
nates before the quantum-classical crossover occurs, so 
the overall resistivity increases with temperature. At 
higher carrier densities, Tp is pushed up to the phonon 
scattering regime, and the quantum-classical cross-over 
physics is overshadowed by phonons so that nonmono- 
tonicity effects are not manifest. 



IV. NUMERICAL RESULTS 
A. Determination of deformation potential 

In this section we provide the numerically calculated 
temperature dependence of the hole resistivity by con- 
sidering both screened acoustic-phonon scattering and 
screened charged impurity scattering. In the calcula- 
tion of phonon scattering we use the parameters cor- 
responding to GaAs: m* = 0.38 me, vi = 5.14 x 10^ 
cm/s, Vt = 3.04 X 10^ cm/s, pm = 5.3 g/cm^ and 
ehi4 = 1.2 X 10^ eV/cm. For the deformation poten- 
tial, we fitted several available mobility data set [Ij, l2^ 
and the best fitted value we obtained is T) = 12.7 eV for 
^depi = (see Fig. [2|). In the following calculations we 
use this value as a deformation potential of p-GaAs. 

To obtain the best fitted value of the deformation po- 
tential, in Fig. [21 we calculate the total mobility p = a/ne 
with Eqs. ([M|) and ([25|) as a function temperature for 
different values of the deformation-potential constant D. 
Knowing the precise value of the deformation-potential 
coupling constant D is very critical because p ~ ^~^, 
i.e., the calculated mobility will be uncertain by a fac- 
tor of four with values of D differing by a factor of 2. 
In this calculation we consider two different scattering 
mechanisms: remote impurity scattering and acoustic 
phonon scattering. We first fit the low temperature data 
(T < 4 K) to find the charged impurity density, nimp, be- 
cause the phonon scattering is severely suppressed and 
the charged impurity scattering determines the carrier 
mobility in this temperature range. We set dimp = for 
simplicity and carried the effect of impurity by nimp- For 
^depi = 0, we find that nimp = 1.22 x 10"^^ cm~^ gives 
the best fitted mobility at low temperatures for the data 
set. With this impurity density we calculate the mobility 
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FIG. 2. (color online). Mobility as a function of temperature 
for several values of D with n = 2 x 10^^ cm~^, nimp = 1-22 x 



10 cm , ridepi = and du 
experimental data |l3 ]. 



0. Black dots represent the 



14 



12-. 



^ 10 



■<^--o.. 



■^■■■o--o., 



0.5 



1.5 



^,,,i(10i2cm-2) 



FIG. 3. Fitted deformation potential D as a function of de- 
pletion density ridepi for a sample with n = 2xl0^^ cm~^ and 
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FIG. 4. (color online). The relaxation rate (solid line) and 
scattering rate (dashed line) as a function of hole energy E 
for (a) T = 10 K and (b) T = 1 K with n = 10^^ cm-^ 
^imp = 0, ndepi = and D — 12.7 eV. DP and PE represent 
the deformation potential and piezoelectric potential contri- 
butions, respectively. 



results qualitatively. 



data at high temperatures (20 K< T <60 K) by chang- 
ing deformation potential. From Fig. [21 we get D = 12.7 
eV as the most suitable value for the p-GaAs acoustic 
phonon deformation coupling constant. 

In Fig. [3] we show the deformation potential coupling 
as a function of the depletion density for a hole density 
n = 2xl0^^ cm~^. The calculated deformation potential 
D strongly depends on the depletion density ndepi for 
^depi < lO"*^^ cm~^, but for higher densities (ndepi > 10"*^^ 
cm~^ it decreases slowly, as seen in Fig. [3l ndepi is a 
measure of fixed charges in the background and typically 
ndepi is unknown. Thus, the uncertainty in the value of 
deformation potential coupling (both for electrons and 
for holes) could be a result of our lack of knowledge about 
^depi- When ndepi = we get D = 12.7 eV. We use 
D = 12.7 eV in the remaining calculations with ndepi = 0. 
However, the different values of D do not change the 



B. Acoustic phonon-limited transport 

Using the theoretical model outlined in Sec. Ill we 
study hole transport limited by acoustic phonons in this 
subsection. In Fig. |4] we show the calculated scatter- 
ing rates r~^ and transport relaxation rates r^^ due 
to acoustic phonons, as a function of the hole energy. 
The relevant transport relaxation rate, r^~^, has been ob- 
tained in Eq. (|26|) . The two characteristic times shown 
in Fig. m differ by the important (1 — cos^) factor |25j . 
The scattering rate r~^ is given by making the replace- 
ment (1 — cos^) ^ 1 in the integrand for the formula for 
r^ given in Eq. (|26|) . rt determines the conductivity (or 
mobility), a = neji = ne^rt/m, where n is the carrier 
density and fi is the mobility, whereas Tg determines the 
quantum level broadening, 7 = fi/2rs, of the momentum 
eigenstates. The scattering time Tg is related to the imag- 
inary part of the single-particle self-energy and simply 
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FIG. 5. (color online), (a) Acoustic phonon- limited resistivity 
of 2DHS as a function of temperature for n = 10^, 10^, 10^°, 
10^^, 10^^ cm~^ and (b) the calculated exponent a in p{T) oc 
T" which is obtained from logarithmic derivatives of (a). 
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FIG. 6. (color online), (a) Acoustic phonon- limited mobility 
of 2DHS as a function of temperature for n = 10^, 10^, 10^°, 
10^^, 10^^ cm~^. (b) Density dependence of coefficient a for 
Tiimp =0 and 5x10^ cm~^, where l//i = l//io + aT in 10 
K< T <60 K range. 



gives the time between scattering events between a hole 
and an acoustic phonon. The difference between rt and Tg 
arises from the subtle effect of the wave vector dependent 
transition rate [25]. The large angle scattering events (or 
large momentum transfer) contribute significantly to the 
transport scattering events, but small angle scattering 
events where cos6> ~ 1 makes a negligible contribution 
to Tt, while all scattering events contribute equally to Tg. 
Our result for the individual DP and total PE rates are 
given in Fig. |4]for n = 10^^ cm~^ at two different tem- 
peratures T = 10 K and T = 1 K. In this calculation 
we take D = 12.7 eV which is the best fitted values of 
the experimental data. We find that rt/rg ~ 1, since the 
screened electron-acoustic mode phonon interactions are 
of relatively short range. It has been known that the 
ratio Tt/rg from remote ionized impurities is much big- 
ger due to the long-range nature of the electron-impurity 
interaction [26|. 

Figure [51(a) shows acoustic phonon-limited resistivity 
of 2DHS in the absence of impurity scattering as a func- 
tion of temperature for different hole densities n = 10^, 
10^, 10^°, 10^\ 10^^ cm-^ The calculated resistivi- 
ties clearly demonstrate the two different regimes: BG 



region characterized by high power law behavior for 
T < Tbg and equipartition region with p ^ T behav- 
ior for T > Tbg- The transition temperature Tbg in- 
creases with density since Tbg ^ V^- As the density in- 
creases the calculated resistivity at a fixed temperature 
decreases. In Fig. [5fb) we show the logarithmic deriva- 
tives of the acoustic phonon limited resistivity, which 
give rise to an approximate temperature exponent of 
acoustic phonon limited resistivity by writing p ^ T^, 
i.e., a = d\np/d\iiT. At low temperature BG regime 
T < Tbg the numerically evaluated exponent a varies 
from 4 to 6 depending on the carrier density. But at 
high temperatures the exponent approaches to 1 as we 
expected, i.e., p(T) ex T. 

Figure [61(a) shows acoustic phonon-limited mobility of 
2DHS in the absence of impurity as a function of tem- 
perature for different hole densities n — 10^, 10^, 10^^, 
10^\ 10^^ cm-^ At high temperatures (T > 10 K) 
the calculated mobilities show very weak density depen- 
dence for the density range n = 10^^ — lO"*^^ cm~^ and 
decrease approximately as p ~ T~^. Thus, the recip- 
rocal mobility increases linearly with temperature, i.e.. 
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FIG. 7. (color online), (a) Resistivity of 2DHS as a function of 
temperature for nimp =0, 1, 2, 3, 5, 7, 10x10^ cm~^ with n = 
10^*^ cm~^, ridepi = and dimp = 0. Dotted lines indicate the 
calculated resistivity due to the charged impurity scattering 
alone, (b) Same as (a) but rescaled by po = p(T = 0.1 K). 



l//i = l/jiio-\-aT, where a is the slope in the relation be- 
tween /i~^ and T. Fig.[6l^b) shows density dependence of 
the slope a for nimp =0 and 5x10^ cm~^ in the temper- 
ature range 10 K< T <60 K. The slope a first increases 
with n, reaches its maximum at n ~ 10^^ cm~^, and 
decreases very slowly for n > 10^^ cm~^. This nonmono- 
tonic behavior is different from that of the n-type GaAs, 
in which the slope a has a minimum value rather than a 
maximum [sl, |27J. 



C. Nonmonotonic resistivity in p-GaAs 

In this subsection we study the transport in the pres- 
ence of both acoustic phonon and impurity scatterings 
and the nonmonotonic behavior in temperature due to 
the competition between these two scatterings. 

In Figs. [3 and [8] we show our calculated total resistiv- 
ity p{T) arising from screened charged impurity scatter- 
ing Pimp(^) and phonon scattering pph(T) as a function 
of temperature. In Fig. [Tl^a) the total resistivity p{T) is 
shown for different impurity densities nimp =0, 1, 2, 3, 
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FIG. 8. (color online), (a) Resistivity of 2DHS as a function 
of temperature for dimp = 0, 5, 10, 15, 20, 25, 30 nm with n — 
10^° cm~^, nimp = 5 X 10^ cm~^ and ridepi = 0. Dotted 
lines indicate the calculated resistivity due to the charged 
impurity scattering alone, (b) Same as (a) but rescaled by 
po-p(T = 0.1 K). 



5, 7, 10x10^ cm-2 with a fixed di^p = 0. In Fig. [H^a) 
the total resistivity p{T) is shown for different impurity 
location from the interface dimp =0, 5, 10, 15, 20, 25, 30 



nm with a fixed impurity density n 



imp 



5 X 10^ 



Figures [ll^b) and E^b) are the same as Figs. [T^a) and 
[HJa), respectively, but rescaled by po = p{T = 0.1 K). 
In a real system the amount of random disorder depends 
on the strength and the spatial distribution of all the 
random impurity scattering centers. However, in these 
calculations we assume that the charged impurities are 
randomly distributed in a 2D plane located at dimp from 
the interface. The calculation is carried out with the hole 
density n = 10^^ cm~^ which corresponds to the Fermi 
temperature Tp ~ 0.7 K. The dotted lines in Figs. [7] and 
[8] indicate the calculated resistivity due to the charged 
impurity scattering alone, Pimp(^)- To calculate the to- 
tal resistivity p{T) we use the total scattering rate of 
Eq. (|M|) because the Matthiessens rule, which is implic- 
itly assumed p{T) = Pimp(^) + Pph(^), is known to be 
not strictly valid at finite temperatures. 

As shown in Fig. [71 when the charged impurity den- 
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sity riimp increases the impurity scattering effects become 
stronger, while the phonon scattering effects are unaf- 
fected. Therefore at high impurity densities the impurity 
scatterings are dominant over phonon scatterings. The 
calculated p{T) increases at lower temperatures (T < 1 
K) due to screening effects, then the quantum-classical 
crossover occurs at the intermediate temperature regime 
around T ~ 1.5 K where nondegeneracy effects make 
resistivity decrease as p ~ T~^. At higher tempera- 
tures (T ^ 10 K) phonon scattering takes over and 
p{T) increases with T. Thus, for large impurity densities 
the temperature dependence of the calculated resistivity 
shows a nontrivial nonmonotonic behavior, arising from 
a competition among three mechanisms discussed above, 
i.e. screening which is particularly important for T < 1 
K, nondegeneracy and the associated quantum-classical 
crossover for T '^ Tp, and phonon scattering effect which 
becomes increasingly important for T > 10 K. At lower 
impurity densities the quantum-classical crossover effects 
are not particularly shown in Fig. [3 because phonon scat- 
tering becomes more important than the classical behav- 
ior p ~ T~^, and the system makes a transition from 
the quantum regime to the phonon scattering dominated 
regime. The linear rise in p{T) for T > 10 K in Fig. [7] is 
the phonon scattering effect. 

The same results shown in Fig. [71 are expected by vary- 
ing the impurity location because the scattering limited 
by the remote impurity becomes weaker as the distance 
of the impurity from the interface increases. In Fig. [8] 
we show the several different kinds of nonmonotonic be- 
havior by varying the impurity location. When the im- 
purities are located very close to the interface (top lines 
in Fig. [8]) the nonmonotonic behavior of the resistivity 
clearly appears in the temperature range we consider (i.e. 
T < 100 K) due to competition among the three mech- 
anisms discussed above. As the separation increases the 
nonmonotonicity becomes weaker because of the reduc- 
tion of the charged impurity scattering and the associated 
weakening of screening effects. In addition, the increase 
of the separation gives rise to the shift of the local max- 
imum peak to the lower temperature. For large separa- 
tions (bottom lines in Fig. [8]) the local maximum peak 
does not appear in the calculated resistivity because it 
shifts to very low temperatures (T < 0.1 K). 

One interesting finding in our calculation is the tem- 
perature region where the calculated resistivity is ap- 
proximately constant, as indicated by the dashed box 
in Figs. [7] and [51 The temperature range of the con- 
stant resistivity appears when the increasing resistivity 
due to phonon scatterings compensates for the decreas- 
ing resistivity due to the nondegeneracy effects. The flat 
region depends critically on the impurity density and the 
location of the impurities, and can be observed in ex- 
periments by varying the doping density and location. 
In Fig. [7] a flat region spanning around 2 K < T < 10 



K appears at an impurity density nin 



3 X 10^ cm-2 



for dimp = 0. In Fig. [8] the flat region for 2 K < T < 
10 K appears at dimp = 5 nm with an impurity den- 



sity nimp = 5 X 10^ cm~^. It is, therefore, possible in 
some situations for a complete accidental cancellation 
between the increasing temperature dependence of the 
phonon-induced resistivity and the decreasing tempera- 
ture dependence of the quantum classical crossover effect 
from impurity scattering in a narrow intermediate tem- 
perature regime. We believe that this has recently been 
observed experimentally ^], but a detailed comparison 
with experiment is not possible due to the complications 
of the parallel magnetic field used in the experimental 
measurement to induce magneto-orbital coupling. 



V. CONCLUSION 

To conclude, we have calculated the temperature 
dependent transport properties of p-type GaAs-based 
2DHSs for temperatures T < 100 K by taking into ac- 
count both hole-phonon and hole-impurity scatterings. 
Our theory includes temperature-dependent screening of 
both charged impurity scattering and phonon scattering. 
We extract the deformation potential D of hole-phonon 
coupling constant by fitting the experimentally available 
mobility data. We find that the deformation potential 
coupling varies (i.e. D = 7.6 — 12.7 eV) depending on the 
value of the depletion density ndepi, which is not known. 
When we assume ndepi = we obtain D = 12.7 eV for the 
p-GaAs acoustic phonon deformation potential, which is 
larger than the generally accepted value in bulk GaAs 
{D = 7 eV) [15] but comparable to the value of the n- 
GaAs (12 - 14 eV) [8, 16, 17] in 2D electron systems. 

We also investigate the nonmonotonicity of p{T) aris- 
ing from the competition among three mechanisms: 
screening, nondegeneracy, and phonon scattering. Both 
screening and phonon scattering mechanisms give rise 
to monotonically increasing p{T) with T (at low tem- 
perature for screening, and at high temperatures for 
phonons), but nondegeneracy effects produce a p{T) de- 
creasing with increasing T for T > Tp. Since phonon 
scattering is the dominant temperature-dependent scat- 
tering mechanism in GaAs holes for T > 5 — 10 K, de- 
pending on the density, the stronger nonmonotonicity 
appears when the impurity scattering is dominant over 
phonon scattering below T ~ 5 — 10 K. We carefully 
study the non-trivial transport properties of p-GaAs at 
the intermediate temperature range (i.e., 2 K< T <10 
K). Interestingly we find that the approximate temper- 
ature independence may appear in which p{T) saturates 
in an intermediate temperature range, arising from the 
approximate cancellation between the quantum-classical 
crossover and phonon scattering. Since this flat region of 
the temperature dependent resistivity depends critically 
on the impurity density and the location of the impu- 
rities, it can be observed in experiments by varying the 
doping density and location. We believe that a recent 
measurement ^] has observed this saturation effect. 
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